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Abstract 

Consider a mixed quantum mechanical state, describing a statistical ensemble in terms of an 
arbitrary density operator p of low purity, trp^ <C 1, and yielding the ensemble averaged expectation 
value tr(p^) for any observable A. Assuming that the given statistical ensemble p is generated 
by randomly sampling pure states IV") according to the corresponding so-called Gaussian adjusted 
projected measure [Goldstein et al., J. Stat. Phys. 125, 1197 (2006)], the expectation value 
is shown to be extremely close to the ensemble average tx[pA) for the overwhelming majority of 
pure states and any experimentally realistic observable A. In particular, such a 'typicality' 
property holds whenever the Hilbert space Ti of the system contains a high dimensional subspace 
Tij^ C 7i with the property that all \ip) G 7^+ are realized with equal probability and all other 
IV") G W are excluded. 

PACS numbers: 05.30.-d, 05.30.Ch, 03.65.-w 



I. INTRODUCTION 



Spheres in high dimensional Euchdean spaces exhibit astonishing geometrical properties, 
as discussed in detail e.g. in basic Statistical Physics lectures: two randomly drawn vectors, 
each connecting the center of the sphere with any point at its surface, are practically orthog- 
onal with extremely high probability; almost the entire volume of the sphere is contained 
within an extremely thin surface layer of the sphere; the latter in turn exhibits an extreme 
concentration of its volume around a very narrow 'equatorial belt', and so on. In quantum 
mechanics, pure states live on unit spheres in Hilbert spaces of usually very high dimension, 
and hence one naturally may wonder about their corresponding peculiarities. One of them is 
the subject of our present paper. Namely, we will show the following main result: Consider 
a mixed state, describing a statistical ensemble in terms of a density operator p with low 
purity, trp^ <^ 1, meaning that the mixed state is very 'far' from resembling any pure state. 
Yet, the statistical ensemble p can be thought of as arising by randomly sampling pure states 
according to some probability distribution. In fact, it is well known (see Sect. II for more 
details), that there are many different probability distributions of pure states which give 
rise to the same mixed state p. Here we show that for any given p of low purity there ex- 
ists at least one such probability distribution with the following quite astonishing property: 
Given an observable A, the expectation value {il)\A\il)) for the overwhelming majority of pure 
states \ip) is extremely close to the ensemble averaged expectation value ii{pA) compared 
to the full range of a priori possible expectation values T[iax.\^){'tp\A\'tp) — m\Ti\^){%lj\A\ip) (the 
latter difference is tacitly assumed to be finite, as is the case for any experimentally realistic 
observable A, see Sect. III). ^ 

For this kind of property, the term 'typicality' has been coined in [1]. While such 'typi- 
cality' results are applicable in principle to general quantum mechanical systems, they are 
obviously of particular interest with respect to the foundation of statistical physics of macro- 



scopic systems at equi 



ibrium, as discussed in detail e.g. in Refs. pj, l2, 13|, IJ, 15|, 16|. Further 



J y, I: 



mm. 



lOj. With our present study we extend previous results from 



related works include 

to yet another important class of probability distributions of the pure states 
I?/'), namely the so-called Gaussian adjusted projected measure (GAP), recently introduced 
in Ref. 
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II. OUTLINE OF THE PROBLEM 



We consider a quantum mechanical system with (separable) Hilbert space Ti of dimension 
N < oo. The system is assumed to be in a mixed state (statistical ensemble) described by 
a density matrix p. Let {\n)}^^^ be an orthonormal basis of eigenvectors of p and p„ the 
corresponding eigenvalues, 

N 

P = ^Pn\n)(n\ (1) 



n=l 



with the usual properties 



Pn > (2) 

N 

J2p- = 1 (3) 

n=l 

In the context of equilibrium statistical mechanics, |n) will usually be the eigenstates of the 
system Hamiltonian, but we will not make use of such a property anywhere in this paper. 
A particularly simple and important example is the microcanonical density operator with 



^Y^w)'<^\ (4) 

or, equivalently, with 



Pmic ^ 



Pn = 1/N+ line S, pn = ^\in^ S . (5) 

where S' is a subset of {1, A^}, consisting of elements (1 < < N). 

Given a density matrix of the general form ([I])-([3]), an arbitrary, normalized pure state 
(e.g. a wave function) can be written in the form 

N 



|^) = ^z„|n) (6) 

n=l 

where Zn '■= {n\ip) are complex coefficients, satisfying the normalization condition 

l|z|| = l (7) 



with the standard definitions 



z := {zi,Z2,...,Zn) (8) 

,n=l / 
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Next, we assume that the statistical ensemble p is generated by randomly sampling pure 
states ([6]) according to some probability density p(z). The corresponding ensemble average 
of an arbitrary function /(z) is denoted by 



/(z):= j dzf{z)p{z), (10) 
where dz represents the natural, uniform measure for the A^-dimensional complex argument 



N 



dz := Y\ d(ReZn) d{lmZn) (11) 



n=l 



For infinite dimensional systems, well defined limits N ^ oo are tacitly taken for granted 
in ( fTOl) and in similar expressions later on. 

Put differently, by averaging pure states (E]), represented as projectors \ip){ip\, according 
to the probability density p(z), the given statistical ensemble p has to be reproduced, i.e. 

N N 

= \n){m\ = \n){n\ = p (12) 

m,n=l n=l 

where the star indicates complex conjugation and where we exploited ([I]) in the last identity. 
Hence, the second moments of the distribution p{z) are fixed by the given statistical ensemble 
p via 

^rn^ri — ^mnPn ■ (13) 

In turn, every p{z) with second moments (II 3p reproduces the given p in (IT^ . We thus recover 
the well known fact that a given density matrix p does not uniquely fix the distribution of 
,.re .a.s Mz). Fo. .no. detaU. expUc. examples, and f.the. .efe.ence. see e., Ref. U. 

Next we consider an arbitrary observable A = : 7i ^ Ti with eigenvectors and 
eigenvalues Ojy, i.e. 

N 

A = ^a,|z/)(z/| . (14) 
i/=i 

Here and in the following we use the convention that Greek labels u and p implicitly refer to 
the eigenvectors of A, which is convenient but somewhat ambiguous in so far as, e.g., |z/ = 3) 
is not the same vector as |n = 3). According to ([1]), the ensemble averaged expectation value 
of A is given by 

N 

(A) := tr(pA) = Y^PnAnn (15) 

n=l 

Amn ■= {m\A\n) . (16) 



Further, any given pure state IV') gives rise to an expectation value (V'l^lV') ^ The 
random distribution of those expectation values, induced by the distribution p{z) of pure 
states, is the quantity of central interest in our present work. In particular, we will be 
interested in identifying conditions under which the expectation values {iIj\A\iIj) will be very 
close to the ensemble average fllSp with very high probability, i.e. for a large majority of 
pure states {ip)- 



III. RELEVANT HILBERT SPACE AND OBSERVABLES 

Without loss of generahty, we assume that the indices n in ([1]) are ordered so that 

Pn > for n < and p„ = for > (17) 

for some integer with 1 < iV+ < A^. The A^+-dimensional sub Hilbert space spanned by 
the basis vectors {\n)}^li is denoted by 7Y+ and the projector onto this subspace by 

N+ 

P+:= . (18) 

n=l 

In particular, P4. is at the same time the identity operator on 7i+ and the microcanonical 
density operator (j4]) takes the form 

Pmic = P+/N+ . (19) 

From ( fT3|) and ( flTl) we conclude 4] that (with probability one) Zn = for n > A^_|_ and 
hence \ip) G H+ according to (jH])- As pointed out in Sect. II, our main goal is to determine 
the probabihty distribution of the expectation values {iIj\A\iIj). It follows that with respect 
to this goal only the restriction/projection 

A+ := P+AP+ (20) 

of the observable A to the subspace is relevant. Equivalently, whenever m or n exceeds 
Nj^ then the matrix element Amn is of no relevance for our purpose and thus can be set to 
zero without loss of generality. 

The full range of possible expectation values {tp\A\tlj) is quantified by 

:= max {iIj\A+\iIj) — min {ilj\A+\ip) 
= max (tplAltp) — min {tplAltp) 

< max(^|A|^) - mm{^\A\^) . (21) 
\tp)en W(^T-i 



The second relation follows from (120|) and the fact that P+ is the projector onto the subspace 
7i+ and the last relation from the fact that 7i+ C Ti. Alternatively, can thus be identified 
with the difference between the largest and the smallest eigenvalues of and is bounded 
from above by the difference between the largest and the smallest eigenvalues of A, cf. Eq. 

Clearly, any given real experimental apparatus has a finite range and hence the corre- 
sponding range of possible expectation values from Eq. fl2T|) is finite. This range is 
specific to the given measurement device, but is (practically) independent of the properties 
(e.g. the size) of the observed system. Here and in the following we restrict ourselves to 
observables A in the above sense. For instance, the energy of a harmonic oscillator is not 
an observable in this sense: in principle, the energy of the oscillator may become arbitrarily 
large (albeit with extremely small probability), but no real device would be able to display 
its value beyond a certain upper limit. Rather, all energies beyond this limit will yield 
one and the same measurement result (e.g. a blow up of the device), and hence only the 
corresponding 'truncated' energy operator would be an admissible 'observable'. 

We emphasize again that the above restriction regarding the admissible observables A 
is of a purely mathematical/formal nature, it does not exclude any relevant observable 
corresponding to a realistic physical measurement. Indeed, it is well known that any realistic 
observable can be built up by means of suitable projector operators, and each such projector 
P only has eigenvalues zero and unity and hence Ap = 1. Essentially, the same 'restriction' 
in fact also applies to the seminal prior works . The measure of distance employed 

in is the trace-norm, involving a maximization over all operators with operator norm 
bounded by unity. Hence, the estimates from [2! become worse and worse, as the maximally 
admissible norm of the considered operators increases. Similar conclusions apply for the 
estimates from The most explicit discussion of this issue is contained in (of. 

To summarize, only the sub Hilbert space H+ and the projected observables ( l20|l are of 
relevance for our purpose. Accordingly, we can and will for simplicity assume in some of the 
following sections temporarily that Ti = 7^+ and thus A = A^. In other words, all subscripts 
'+' will be omitted and it will be taken for granted that p„ > for all n. Furthermore, we 
can and will focus on observables of finite range according to (1211) . 
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IV. THE GAUSSIAN ADJUSTED PROJECTED MEASURE (GAP) 



To avoid unnecessary technical complications, we temporarily restrict ourselves to finite 
dimensions A^. In the final results of our calculations, the limit N ^ oo can be readily 
performed. Furthermore, we assume p„ > for all n without loss of generality, see end of 
Sect. III. 

Taking for granted the above assumptions that N < oo and p„ > 0, we define 

p{z) := / rfyATexpl-f;^) | |y| p 5(z - y/| |y| |) . (22) 

V n=l ^" / 

Here, dy is defined like in (fTTj) . i.e. we are dealing with an integral over 2N real integration 
variables, and Af := Y[n=i('^Pn)~^ is a normalization constant (see below). The delta- 
function is by definition understood in such a way that the ensemble average of an arbitrary 
function /(z) from flTU]) takes the form 

m = [ dyAfexp f-f;^ ) ||y|r/(y/||y||) (23) 

V „=i J 



Eq. ( l22i) defines the Gaussian adjusted projected (GAP) measure [4!] associated with the 
density matrix ([1]), written in the form of a probability density p{z) with respect to the 
natural measure (ITT!) . The word 'Gaussian' in the notion GAP refers |J] to the exponential 
factor in fl22|) : the word 'adjusted' refers to the factor | |y | P, which is needed to fulfill condition 
( IT2l) (see below); the word 'projected' refers to the delta- function in ( !22l) . guaranteeing the 
normalization condition ([7]) (see below). 

In the special case of a microcanonical density operator we have pn = ^/N for all n 
according to ([5]) and our assumption at the beginning of this section, yielding with (123|1 the 
result 



/(z) = /^/yAre-llyll'^||y|r/(y/||y||). (24) 



"or arbitrary unitary N x N matrices U. Since also /(z) is 
J] that two arguments z with equal length ||z|| are realized 



It follows that f{Uz) = /(z) 
arbitrary, we recover the fact 
with equal probability. 

Returning to the general case, it is often convenient to change from a Cartesian represen- 
tation of the complex numbers y„ in terms of real and imaginary parts (cf. ffTTj) ) to a polar 
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representation in terms of r„ > and (pn G [0, 2ti) via the usual relation ?/„ = r^e*'''". Then, 
the ensemble average of an arbitrary function /(z) from fl23|) can be rewritten as 

TV 



/(z) 



n 



n=l 



oo /•27r ^ 

Jo T^Pn 



7(c) 



(25) 
(26) 



We first consider the special choice /(z) := 1, implying with fllOj) and (123]) that 



dzp{z) 



■ TV 

n 



oo /•27r 



N 
n=l 



(27) 



The integrals over ipi are trivial, each yielding a factor 27r. Hence, we can infer that 

N N 



* ^ * /"OO 

dzp(z) = x^n / 



71=1 1 = 1 

The integrals over for / ^ n are of the form 

2ri 



Pi 



(28) 



dri 



V «n. 



1 . 



(29) 



Likewise, the integral over ri for Z = is of the form 



Pi 



d \ _„2 / 



Pi 



(30) 



where the second identity follows by a partial integration and the last identity by means of 
(^^. All in all, the right hand side of (|27|) thus amounts to J2n=iPn with we see 
that p{z) is normalized to unity. Observing that the right hand side in (122|) is non-negative 
for any z, we can conclude that p(z) is indeed a well-defined probability density. 

Next, we consider the special choice /(z) := 6 (^X — J2n=i k^P) arbitrary real 



number X. According to fl2B]) . the argument /(c) in fl2^ takes the form 6{X — 1) and 
thus can be brought in front of all the integrals. The remaining integral is identical to 



the one evaluated in the preceding paragraph, i.e. its value is unity, and hence /(z) = 
6{X — 1). It follows that p{z) indeed takes non-zero values only for arguments z respecting 
the normalization condition ([7]). 

Finally, we consider the special choice /(z) := z^Zn. Exploiting fl25|) we obtain 



ri 



dri / d(pi — e 
io T^Pi 



n 

i=\ 



m.' n.^ 



(31) 
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If m 7^ the integral over can be carried out first, being proportional to Jg ^ dip^ e*"^" = 0. 
Hence z'^Zn = if m 7^ n. In the case m = nwe have e*'^"'^™"'"'^") = 1 and we are left with 
independent integrals of the form J^^ dip = 2ti. The remaining integrals over r; are of the 
same type as those already encountered in ( l29l) and (!3Ql) yielding the final result z'^Zn = Pn- 
All together we thus find that the GAP measure ( l22l) indeed fulfills the condition ( |T3l) and 
hence reproduces the correct statistical ensemble (fT2l) encoded by the preset density operator 
p from ([1]). Without the 'adjusting factor' | |r| p in (!22l) this property could not be maintained 

3- 

V. EVALUATION OF THE VARIANCE 

As pointed out in Sect. II, our main goal is to determine the probability distribution of 
the expectation values {iplAl^p) induced by the distribution p(z) of pure states according to 
the GAP measure (|22l) . For the first moment, {iplAl^p), the expected result {A) is readily 
recovered by means of f[T21) and f[T^ : 

N 

(VPM = tr(|^)(^|A) = triWMA) = trip A) = Y,Pn A^n = {A) . (32) 

n=l 

In the present Section, our focus is on the variance 

4 := mA\i^) - {A)f = imi^)' - ¥\m' ■ m 

We emphasize, that this variance characterizes the dispersion of the expectation value of A 
for different pure states {ip), and not the "quantum fluctuations" associated with individual 
measurements of A of a fixed pure state {ip). 

Observing that the two observables A and A — (A) have the same variance and the same 
range according to ( 12T1) . we can and will restrict ourselves in this section without loss of 
generality to observables A with the property 

N 

Y,PnAnn = {A) = . (34) 

n=l 

Furthermore, we maintain the assumptions N < 00 and p„ > for all n, as intro- 
duced at the beginning of the previous Section. In particular, we thus have A = A^ 
and both the eigenvalues (cf. dHD) and the diagonal matrix elements Ann (cf. ( TT6i) ) 
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are bounded from above by amax '■= niax|^)g7^('?/'|y4|'?/') = max^ai, and from below by 
dmin '■= ^^^\i})eH{4'\^\4') = minyajy. In view of fl34l) it follows that amax > and amin < 
and hence with (^11) that 

I -^n/n I ^ I f^iy I 



^nnh I 



With of (ED and flMI). the variance fl33l) takes the form 



AT TV 



j,k=l m,n=l 



(35) 



(36) 



The average in the last term can be rewritten by means of (!25|) as 



n 

.1=1 



2n 

Trpi 



dri I dipi — e-^'/P' 



J' k' m' 



i{-lpj+ipk-'Pm+'Pn) 



(37) 



The evaluation of these integrals is analogous but somewhat more involved than those from 
the preceding Section: The integrals over the angles ipi can be readily performed, yielding a 
factor of (27r)^ in the two cases (i) j = k and m = n, (ii) j = n and k = m, and zero in any 
other case. Taking care not to count the case j = k = m = n twice and after a convenient 
renaming of the summation indices we thus obtain 



fr r dn ^ e-^f/^' 

In order to evaluate the integral Imn, we consider the auxiliary function 



^2 ^2 

m n 



(38) 
(39) 



/i(x,y) 



■ N 

n 

.«=i 



2r, 

Pi 



N 



exp <^ - ^(x + yz)? 



(40) 



1=1 



where x > 0, y := {yi, ....,yj\f), and > for all n. Observing that the right hand 
side in ( 1401) factorizes into independent integrals of the form dririe~^^ = 1/26 with 



h := X + yi > Q (see also fl29l) ). we obtain 



f^^ Pi x + yi 



(41) 



Next, we note that the integral over the x-dependent terms in f l40|) is of the form 
(ix e~^ll'"ll^ = l/||r|p (see (E])), implying that 

N 



H{y) := / dxh{x,y) 
'o 



n 

.1=1 



, 2ri 
dri — 

Pi 



N 



exp< - ^yirf W||r| 



(42) 



1=1 
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By comparison with (l39l) we can conclude that 



'IfL'n r\ o 

oymoyn 

By combining (HT|) - (l43|) it follows that 



(43) 



I = I dr 



oo q2 ^ 







n-— 

11 T), -r -I- 



= PmPn{l + Smn)Kmn (44) 

?/i=l/Pi 



dymdyn fJ^Pi x + yi 

POO 

Kmn ■= / dxgmn{x)G{x) (45) 

Jo 

gnm{x) ■■= {1 + Xp^y^ {1 + XPny^ (46) 
N 

G{x) := n(l + ^P')"' (47) 
1=1 

Finally, this yields for the variance ( l38i) the result 

AT 

(^A ^ ^ {^mm-^nn ~l~ ■^mn^nmlPmPn-^mn ■ (48) 

m,n=l 

Next we turn to a more detailed discussion of the integrals Kmn in (H5l) . Clearly, the 
integrand is a positive function of x, bounded from above by unity, and decaying like 
for large x due to our assumption that p„ > for all n, see below (IMl) . Hence the integrals 
Kmn are finite and positive. Specifically, for the microcanonical density operator we have 
Pn = for all n according to ([5]) and our assumption below yielding with ( H5l) the 
exact result 

Kmn = I dxil + x/Ny''-' = . (49) 







To further eva 
Taylor's theorem 



uate Kmn in the general case, we rewrite gmn{,x) from ( H6l) by means of 



12| as 



Qmnix) = gmn{0) + Xg'mniO) + —gmn{x9mn{x)) 

= 1-X{pm+Pn) + X^ {Pm + PmPn + pI) Xmn{x) (50) 

for certain functions Omn{x) and Xmn{x), satisfying 9mn{x), Xmn{x) G [0, 1] for all x > 0, but 
for the rest depending in a non-trivial manner on x, m, and n. Note that while an infinite 
power series expansion would not converge for arbitrary x > 0, the above finite order Taylor 
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expansion is an exact identity [12] for all x > 0. As a consequence, (H5l) can be rewritten as 

= K^^) -(p^+p^)K^^) +2{pl+p^p„,+pl)K,nnK^^^ (51) 

K^^) ■= / dxx''G{x), A; = 0,1, 2 (52) 
k' Jo 

Kmn e [0,1] (53) 

From (1471) and p„ > for all n we can infer that the integrals in flS^ are finite (and positive) 
if and only if 

N>4. (54) 

The latter condition is tacitly taken for granted henceforth. 

Next we return to the variance in (HHll . In view of we see that sums of the form 

N N N 

^ ^ -^mm^nnPmPnQmn ^ ^ -^mrriPrn ^ ^ -^nnPnQ mn (55) 
m,n=l m=l n=l 

are zero if the coefficients Qmn are either independent of m or independent of n. Hence the 
first sum on the right hand side of fHS]) vanishes in the special case corresponding to the 
microcanonical density operator. Likewise, in the general case we can conclude with ( 15T]) 
that the first sum on the right hand side of (l48l) takes the form 

N N 

^ ^ -^mm-^nnPmPnK-mn 2,K^ ^ ^ ^ ^mm^nriPmPn (Pm, ~^ PmPn ~^ Pn) ^mn ; (56) 
m,n=l m,n=l 

yielding with (1531) and fC^^-* > the estimate 

N N 

I ^ ^ ^mm^nnPmPnKmnl < 2i^'^^^ ^ ] I ^mm | | ^nn | (2pmP^ + PmPn) • (57) 
m,n=l m,n=l 

With (ED and (ESI) we obtain 



N I ^ 

^ ^ ^mm^nnPmPnKmn ^ 27^^^ ''^yl j ^ ^ ^ Pn ~^ 



N 
n=l 



(58) 



m,n=l \ n=l 

Turning to the second sum on the right hand side of (j48|) . we note that 

N N 

^ ^ ^ ^mn-^nmPmPn^mn ^ ^ ^ ^ ■^mn-^nmPmPn • (59) 

m,n=l m,n=l 

The first inequality follows from the fact that Amn^nm = l^mnP > 0, PmPn > 0, and 
Kmn > for all m, n. The second inequality follows from gmn{x) < 1 according to (H6l) . 
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hence Kmn < K^^^ according to fHS]) and flS^ . By means of ([T]) and ffTBl) one readily finds 
that 

N N N N 

tr(py4)^ = ^(m|py4^ \n){n\pA\m) = ^ Pm{m\A\n)pn{n\A\m) = ^ PmPnAmnAnm ■ 
m=l ?i=l m,ri=l m,n=l 

(60) 

Likewise, by using the eigenvectors \i>) and eigenvalues a^, of A from (fT4|) to evaluate the 
trace one obtains 

TV Af N N 

tr{pAf = '^{fi\pA'^\u){u\pA\fi) = {p\p\u)a^{u\p\p)a^ = ^ a^^a^p^^p^^ . (61) 

(1=1 l/ = l IJ,,V=1 fl,U=l 

Combining ([59])- ([M]) and Pf,uPu^l = IPfiA"^ yields 

N N N 

J2 A^rAmnPmPnKmn < J] a,a^\p^,\'' < K^"^^ ^ \p^,\' = iT^o^Ajtrp^ , (62) 

m,n=l fi,u=l fi,iy=l 

where the second inequality follows from fl35|) and the last equality from fl6T]) with A = 1. 

In the special of a microcanonical density operator (j4]) we have seen below (!55|) that the 
first sum on the right hand side of ( HHI) vanishes. Exploiting ( H9l) . ( l60l) . and the fact that 
p„ = 1/A^ for all n we obtain 

N 



N , ... trA2 .fr^ 



= , -I tr(pmicA)" 



(63) 



AT+i ^ A^(A^ + 1) Ar2 

where the eigenvalues of A, see (IT^ . 

Returning to the general case, the variance fHHl) can be estimated from above by means 
of (158|) . fl62|) . and the relations ^Y^.=\Vn — ^'^P^i J2n=iPn — (trp^)^/^, derived in Appendix 
A, as follows 

< K(°)Altrp2 + 2K^^^Al (2 [tip^f + [trp^]^) . (64) 

Our next goal is to find upper and lower bounds for G{x) for x > (a; = is trivial) in 
order to estimate K^''^ from (1521) . To this end, we consider x > as arbitrary but fixed, and 
consider the right hand side in flTTI) as a function of p := (pi, ...,pn), 

N 

g(p) := Hil+xpn)-' . (65) 

n=l 
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The basic idea is to determine its maximum and the minimum under the three constraints 
([2]), ([3]), and pn < Pmax for all ra, where 

Prnax ■= maxp„ . (66) 

n 

The differential/variation of ( l65l) reads 

6Q{p) := -a;g(p)f;-4^, (67) 

complemented by the constraints J2 ^Vn = 0, 6pn > if Pn = 0, and 6pn < if p„ = Pmax- 
Observing that xQ{p) > on the right hand side of fl67|) and that the factors 1/(1 + xpn) 
are smaller (but still positive) for large p„ than for small p„ implies that Q{p) can always 
be made larger {SQ{p) > 0) by making the already large pn still larger {6pn > 0) and the 
already small Pn still smaller {6pn < 0). As a consequence, Q{p) is minimal if all p„ are 
equal, implying that 

G{x)>{l + x/Ny^ . (68) 

On the other hand, Q{p) cannot be increased any more if and only if the small p„ have 
reached the lower limit pn = and the large p„ the upper limit pn = Pmax- Denoting by 
Nmax the number of those p„ equal to Pmax, their total weight N^axPrnax is generically still 
not exactly equal to unity for any integer N^ax- Hence there must remain one weight Pn 
with a value 1 — NmaxPmax ='■ Po £ [0,Pmax] ordcr to fulfill the constraint ([3]). All in all, 
this implies the upper bound 

G{x) < (1 + xp„„,)-^— (1 + xpo)-' . (69) 

Next we note that for any a > the auxiliary function f{y) := ln(l + ay) — yln{l + a) is zero 
for y = and y = I and has a negative second derivative for y > 0, implying that f{y) > for 
all y e [0, 1]. Setting a = xpmax and y = Po/pmax it follows that (1 + xpo)^™" > i'^ + xprnaxY" 
and due to po := I - N^axPmax that 1 + xpo > (1 + xp maxY ^ ~ ■ With ([69]) we thus 
can infer that 

G{x) < (1 + a;p„„,)-^/P— . (70) 
From ([3]) and (!66|) we see that 

Pmax > l/N . (71) 
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Further, the upper bound (1701) yields finite integrals (152|) only if 



Pmax < 1/3 . (72) 

Note that this condition implies > 3 and hence condition (l54l) is automatically satisfied. 
Taking for granted fl72|) inequality we can infer by exploiting the bounds fl68l) and flTOj) in 
(!52|) and after performing A; partial integrations that 



^^'^ = rilT^' ^ = 0'1'2 (73) 
G [l/iV,p^,,] (74) 



With fl64|) we thus obtain for the variance the upper bound 

. < A2 ^ trp^ 4 [trp^/^ + 2[trpT \ ...^ 

\ fmax V l/max )\^ ^l/max j V ^l/max ) / 

VI. DISCUSSION OF THE MAIN RESULTS 

The upper bound ( l75ll for the variance from ( l33ll is the first main result of our paper. 
In our derivation we have assumed that {A) = (see but since the variance from 

fl55]) and also all the other quantities appearing in fl75]) remain unchanged upon replacing 
Ahj A — (A) we can conclude that fl75|) remains valid for arbitrary A. Moreover, we made 
the assumption that N < oo and p„ > for all n in deriving (!75|) . Since neither of the 
quantities appearing in the final result fl75l) give rise to any problem in the limit p„ —>■ 0, 
the assumptions p„ > can be given up as well. Finally, the limit ^ oo depends on 
the meaning and existence of this limit for the quantities appearing on the right hand side 
of (|75]) . In particular after dividing both sides by (see below), we expect that in many 
important cases this limit will not give rise to any problems. The only remaining condition 
for fl75|) to be apphcable is thus Pmax < 1/3 (see fl72l) ). 

In the special case of a microcanonical density operator (j4]) we have obtained as a second 
main result the exact relation ( 163|) for the variance under the same assumptions as above, 
namely (A) = 0, N < oo, and p„ > for all n. Accordingly, for more general observables 
A with (A) 7^ we have to replace A hj A — (A) in (jHS])- Next, if not all p„ are positive 
and thus equal to 1/A^, then we have to replace A by A+ and A^ by A^+, as discussed at the 
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beginning of Sect. III. All in all, we thus obtain for a microcanonical density operator 
the general exact result 



2_ +^r. (A I A \M2 _ I "^^4 



o 1 / trA 



1^ 



2^ 



(76) 



where is the projection of the original operator A onto the subspace spanned by the 
basis vectors \n) with non-trivial weights p„ > 0, see Eq. (!20|) . As before, the meaning and 
existence of the limit N —>■ oo depends on the behavior of A^+, try4_|./A^+, and trA'^/N^ in 
this limit, but is expected not to give rise to any problems in many important cases. 

Results similar to ( 1761) have been previously derived in Ref. in Ref. ^ (see formula 
(C.17) therein), and in Ref. [ll| (see Lemma 3 therein, whose proof is very close in spirit 
to Ref. [l|]). The main difference is that these results only apply to the special case that 
5 = {1, .■.,N} and hence = in ([5]), implying that pmic in & is proportional to the 
identity operator, cf. (fT9l) . At first glance, a further difference appears to be that the above 
mentioned results do not refer to the GAP measure ( l22l) associated with the above pmic 
but rather are derived under the assumption that all (normalized) pure states & Ti. are 
realized with equal probability. However, by noting that the latter assumption uniquely 
determines the probability density p{z) for the coefficients z„ in ([6]) and that the GAP 
measure does fulfill the assumption (see below Eq. (!24l) ) we can conclude that there is in 
fact no difference in this respect. As a by product we can infer that ( l76l) in particular applies 
to the case that all (normalized) pure states {ip) within the subspace 7i+ C 7i are realized 
with equal probability and all other \ip) & Ti. are excluded. After submission of this paper, 
A. Sugita pointed out that the same finding is also contained in his recent work j^. 

Of particular interest in (1751) are situations for which the bracket on the right hand side 
becomes a small quantity. Therefore, we now focus on the case that the so-called purity trp^ 
of the mixed state p is low, i.e. 

trp2 < 1 . (77) 

We recall the well-known facts that the purity is one if and only if p corresponds to a pure 
state (p = for some {ip) G 7i), is smaller than unity in any other case, and takes the 

minimal possible value 1/ N if Pn = for all n in ([1]). Roughly speaking, a low purity 
trp^ thus means that the mixed state p is very 'far' from any pure state {ip) G Ti. 

According to Appendix A, the quantity Pmax from (!66l) can be estimated from above and 
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from below as follows 

trp^ < pmax < A/trp2 . (78) 

Hence, assumption fl77j) is fulfilled if and only if all p„ in ([1]) are small and is tantamount to 
the condition 

Pmax -C 1 . (79) 

Accordingly, there cannot be just a few dominating p„ in ([T]) in the sense that their sum 
would already be of the order of unity. In particular, the dimensionality of the Hilbert 
space H must be large according to ( 17T1) . 

Exploiting (1771) and fl78|) in fl75|) yields our final main result 

a2 < Altrp2 (1 + 0{^/t^)) . (80) 

By rewriting fl76l) in a form analogous to the last expression in fl63|) one readily sees that 
the upper bound (IHOjl is expected to be quite tight in typical cases. 



VII. SUMMARY AND CONCLUSIONS 



Given any mixed state p of low purity ( 1771) there exists at least one probability distribu- 
tion p{z) of pure states ([6]), namely the GAP measure (l22l) . with the following properties: 
(i) By randomly sampling pure states \ip) according to this probability distribution, the 
preset statistical ensemble p is reproduced, (ii) Given any observable A, for the overwhelm- 
ing majority of pure states {ip) sampled according to p{z) the expectation value ('?/'|y4|'?/') 
deviates extremely little from the ensemble averaged expectation value tr(pA) compared 
to the full range of a priori possible outcomes of a measurement corresponding to A. 
The latter statement can be expressed more rigorously by means of (IHOj) in combination 



with Chebyshev's inequality [l3|], and it is tacitly assumed that this range Aa from fl2T|) is 
non-zero and remains bounded even in the case of an infinite dimensional Hilbert space Ti., 
as is the case for any experimentally realistic observable A (see Sect. HI). 

On the one hand, in general there are other measures p{z) besides the GAP measure 
which also satisfy property (i) above but for which property (ii) may not necessarily remain 
true. On the other hand, p fixes all observable properties of the system via ( JT5l) . so that 
under typical circumstances any further information regarding p{z) is neither necessary nor 
available. Hence, in order to uniquely specify p(z) for a given p, one has either to introduce 
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and justify additional postulates regarding p{z) , or to show that many or all of the 
p{z) compatible with p lead to essentially the same final conclusions (ii) |5], or one has to 
include the preparation and equilibration process of the system into the consideration j^. 
In our present work we have focused on the first among those three options. ^ 

The justification for selecting the GAP measure has been discussed in detail in Ref. [^]. In 
particular, it is argued in that this measure arises naturally when considering macroscopic 
systems in thermal equilibrium and hence is the most appropriate choice, at least in cases 
when p is known to be the canonical density matrix. Furthermore, as shown in Ref. {4] and 
again in Sect. VI, this measure is the unique solution in the case of a microcanonical density 
operator (jlj) under the additional assumption that all (normalized) pure states {ip) G 7^+ 
are equally likely and all other \ip) E H are excluded, where 7i+ is the sub Hilbert space 
spanned by all the eigenvectors \n) with pn > in (Q, i.e. the quantum mechanical analogue 
of the classical energy shell within the standard microcanonical formalism. 

In other words, whenever the Hilbert space 7i of the system contains a subspace C 7i 
with the property that all {ip) G 7i+ are realized with equal probability and all other G Ti 
are excluded then the variance A^, characterizing the dispersion of the random variable 
{ip\A\ilj) (see fl33|) ). is given by the exact relation fl76|) . where is the dimension of 
and the restriction/projection of A to 7Y+ (see ( 120|) ). 
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APPENDIX A 

By means of ([1]), the definition of Pmax in f l66|) and the normalization ([3]) we can conclude 
that 

N N 
n=l 71=1 

for any integer k > 0. For the microcanonical density operator (jl]), the above inequality 
becomes an equality, i.e. the lower bound for Pmax following from (1A1I) cannot be improved 
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in general. Likewise, one readily sees that 

TV 

trp*^ = > maxp^ = v^a. ■ (A2) 

n=l 

Here the inequality becomes an equality if — > 1 for one index n and for all 

m ^ and hence again no general improvement of the corresponding upper bound for Pmax 
is possible. In particular, for /c = 2 we have p^ax — ^rp^ and hence we can conclude that 

< < pL. = {Plaxf' < itrp'r^' (A3) 
for any integer k > 1. Finally, this result yields 

N N 

trp' = Y.Pn< J2p'm-alpl = pL>p' < {i. p't-'^ 'h. = {i.p'f' (A4) 
n=l n=l 

for any integer k >2. 
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